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FOREWORD 

This report presents the results of a study of 
SOM of the therMdynaaic and fluid Mchanics aspects 
pertaining to the aanageMnt of cryogenic fluids in a 
systen in which settling of the liquid and gaseous 
phases of the fluid is achieved through the use of a 
dynauic device (a helical screw) located within the 
propellant tank. 


The study was perfomed under contract HAS9- 16549. 


ABSTRACT 


The settling behavior of the liquid and gaseous 
phases of a fluid in a propellant tank in a saro-g 
environnent, when such setf.ling is induced throu|^ tha 
use of a aynanic device, in this particular case, a 
helical screw was studied. Particular Mphasis was 
given to: (1) the description of a fluid nachanics 

nodel which seems applicable to the system under con- 
sideration, (2) a First Law of Thermodynamics analysis 
of the system, and (3) a discussion of applicable scaling 
rules. 


a 


a 


a 


iv 


TABLE OP CONTENTS 


Ptft 

Poreword . ill 

Abstract iv 

Table of Contents ..... v 

List of Pigures yi 

Introduction 1 

Technical Discussion 

A. General Considerations 4 

B. Development of a Model 

Suitable for the Description .... 

of t^ Fluid Mechanics 

Behavior of the System 7 

C. Thermodynamic Considerations .... 23 

D. Power Requirements 47 

E. Comments and Computational 

Procedure 56 

F. Similarity Analysif and Model . . . 

Testing 62 

Conclusion 71 

References 72 


V 


LIST OF FIGURES 


Page 


Figure A-1 

Helical screw- tank 
arrangeBent 

4 

Figure A- 2 

Helical screw 

S 

Figure B>3 

Relationship between 

^f/t ^s/t 

9 

Figure B-4 

Velocity diagraa for the 

case in which - 0 .... 

10 

Figure B*S 

Velocity diagran for the 
case in which v . vl .... 

^f/t ^ 

11 

Figure B-6 

Velocity diagram for the 

case in which V > V^ 

12 

Figure B-7 

Velocity diagram for the 

case in which 1 > S > 0 .... 

IS 

Figure B-8 

Velocity diagrams for 
particles at r^^, r£ and r^ . . . 

14 

Figure B-9 

Velocity diagrams at r ■ 0 

and r • R. 

s 

15 

Figure B-IO 

Distribution of the tangential 
velocity components of the 
fluid 

17 


Vi 



List of Figures (Cont.) 


Ptf# 


Figure B-11 

Velocity components for e 
particle of fluid in 
cylindrical coordinates 
(r, e, i) 

18 

Figure B-12 

Settling configuration . . . 

21 

Figure C>13 

Control volume 

23 

Figure C-14 

Fluid configuration for 
rotational kinetic energy 
considerations 

39 

Figure C-15 

Differential element of 
mass 

42 

Figure D-16 

Flow configuration between 
the tip of the screw and 
the tank wall 

49 

Figure D-17 

Channel -like conduit through 
the helical screw 

53 

Figure D-18 

Relation between P, a, and s. 

54 

Figure E-19 

Inside of the screw shaft 
used as a return path for 
the liquid 

57 


vii 


I 


INTRODUCTION 


The advent of space flight and the problems idilch 
have arisen in the management of cryogenic fluids in low 
or other adverse g conditions, in particular, those 
concerning withdrawal of the proper phase of a fluid, 
either the liquid or the gaseous one, from a propellant 
tank, have led to the proposal of several solutions to 
the problem. 

Some of the proposed ideas have led to the develop- 
ment of hardware which has demonstrated the capabilities 
and feasibility, as well as the limitations of the 
concepts, especially in view of the application for idiicu 
they were intended. 

Engine restart, for example, requires some sort of 
device to retain and subsequently supply enough liquid 
to the engine's feed line during the initial moments idien 
the engine is fired, until the main body of liquid is 
settled by the thrust force. 
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However, it is probably undesirable to rely on 
thrust forces to achieve settling of the phases for the 
purpose of transfering large quantities of cryogens 
between tanks in space, because of the tine required by 
the transfering process. 

Passive devices (those needing no external neans 
to achieve the outflow of the required phase) , especially 
of the screen type, will probably be the first choice for 
a transfering process. Certain devices of the dynanic 
type (those driven by an external energy source) which 
have been proposed, look as attractive alternatives for 
this purpose. 

This study analyzes the behavior of a systen in 
which the settling of the phases is achieved through the 
use of an externally driven device located inside the 
propellant tank, and establishes procedures for the 
determination of some of the more intevesting parameters 
of the system. The device used is a helical screw. 

The tank is cylindrical with hemispherical ends, and the 
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«xis of rotation of the screw lies along the longitudinal 
axis of the tank. The screw runs from one end of the 
tank to the other with a very little clearance between 
the tip of the screw blade and the tank's wall. 

The external energy supplied to the screw is 
transferred to the fluid through body forces. Soue of 
the energy is dissipated through frictional effects due 
to the viscosity of the fluid. 
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TECHNICAL DISCUSSION 


A. GENERAL CONSIDERATIONS 


The system studied, consists of a helical screw 
rotating within a cylindrical tank (Figure A-1). Power 



Fig. A-1. Helical screw- tank arrangement 


to rotate the screw is supplied from an external source. 
The ends of the tank are hemispherical. 

The length of the cylindrical portion of the tank 
is denoted by while the radius of the hemispherical 
ends is denoted by R^. The latter is of the same length 
as the radius of the cylindrical section. 
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The following parameters characterize the helical 
screw (Figure A-2): diameter, helix angle, pitch, flight 

nimber, flight depth and screw length. 






Fig. A-2. Helical screw 

The screw diameter, D^, is the diameter of an 
imaginary cylinder into which the screw would closely fit. 
The angle at which the helix cuts the normal circular 
section of this cylinder is the helix angle a. A flight 
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is one conplete turn (360^) of the screw helix. The 
pit:h, P, is the linear distance between adjacent flights. 
The screw length is determined by the longitudinal length 
of the tank since the helical screw should extend from 
one end of the tank to the other. The flight depth, fd, 
is the radial distance measured from the root to the tip 
of the helix. The helix is a right one which means that 
a line, such as A-B, drawn in a radial direction and 
lying on the surface of the screw, is perpendicular to 
the longitudinal axis of the screw. (Helical screws in 
which the line A-B makes an angle with respect to the 
axis different from 90** are also possible, but only a right 
one will be considered in the analysis which follows) . 


# 


B. DBVBLwPNBNT OF A NODBL SUITABLB FOR THB DBSCRIPTION 
OF THE FLUID MBChANICS BBHAVIOR OF IHB SYSTEM. 

As the helicml screw turns at the rate of N 
revolutions in unit tine, a point in a fixed "nut" 
engaged to the screw would advance a distance equal to 
the pitch F in each revolution of the screw. Thus, the 
rate of advance V* of the point would be 

V‘ - NP (Bl) 

However, a fluid particle would advance in a 
direction parallel to the axis of rotation of the helical 
screw with a speed V which would depend upon the slip, 
defined as the excess of V* over V. Defining the ratio 
of the slip to V as the slip ratio, S, then 

s - VL-_V . 1 . Y . >1, (B2) 

V V NP 
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i. e. , 


Introducing the pitch to dieneter ratio, X, 



(BS) 


then, 

S . 1- X— (B«) 

NXDs 


Defining the advance ratio, J, as 


J - 



then, 


S - 


1 - 


J 

X 


(BS) 


(B6) 


J, S, and X are non*dinensional quantities and fora 
part of the scaling laws of the system. 

The velocity of a particle of the fluid can be 
described by either a frame of reference fixed with res* 
pect to the tank or by a frame of reference fixed with 
respect to the screw. 
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Let be the velocity of a particle of the 

fluid as seen by an observer located on the rotating 
screw. Since adjacent flights of the screw provide a 
channel through which the fluid can flow, the velocity 
of the fluid, would be directed essentially along 

lines parallel to the surface of the screw, as shown in 
Figure B-3. 


Let be the velocity of the particle as seen 

by an observer fixed with respect to the tank, while 
Vg/£ is the velocity of a point in the helical screw as 

seen by the same observer. The latter velocity is 
always tangent to the circle traced by the point as the 
shaft rotates. 



axis of rotation of 
the helical screw 


I 



Fig. B-3. Relationship between V£^^ 


''{/t ^/t 
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The relationship between ^f/t* ^s/t 

is shown in Figure B-3 and is given by 


'{/t 


\/t * Vf/, 


(B7) 


Several flow conditions will be examined as follows: 


a) Suppose • 0- The corresponding velocity 

diagram is shown in Figure B-4. In this case, 
the particle does not advance in the direction 
of the axis of the screw, but performs a rigid 
body rotation with the screw. Thus V ■ 0 and 
S - 1. 

r oxis of rofotion of 
the helicol screw 


f 


=» Vf/, • V»/( 


Fig. B-4. Velocity diagram for the 
case in which 


Consider the situation in which « V. 

The corresponding velocity diagrta is Shown 
in Figure B-5. In this case V - V* and S ■ 0. 
The particle's advance along a direction 
parallel to the axis of the screw corresponds 
to the advance that a point in a fixed 'nut' 
engaged in the screw would experience. 



Fig. B-5. Velocity diagrga for^the 
case in which Vf/^" V. 
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c) Figure B-6 shows the velocity diegru for the 
case in which V > V* so that S < 0. 



Fig. B-6. Velocity diagram for 

the case in which V > V 


The magnitude of and of &re larger than 
those of the previous case for the same value of V. 

d) The velocity diagram, for the case in which 

and are smaller than those in (b) for the 

same value of V, is shown in Figure B-7. The 
value of S lies between 1 and 0, approaching the 
latter value as the particle approaches rigid 
body rotation about the axis of the screw. 
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QKis of rototion of 
tht htlicol tcrtw 


I 



Fig. B-7. Velocity diagram for the 
case in which 1 > S > 0 

Thus far, the velocity diagrans for a single particle 
have been considered. We will now consider the velocity 
diagrams for a set of particles lying on a radial direc- 
tion. Since the tangential velocity of points in the 
screw relative to the tank, varies linearly with 

radius, it is inmediately apparent that we will have 
different velocity diagrams at different radii. Further- 
more, these velocity diagrams will depend on the magnitude 
of the axial component of the velocity of the flu.d with 
respect to the tank (i. o. V) for all points of a cross 
section normal to the longitudinal axis of the tank. 


QhiGi'iAL IS 
OF POOR gUALOY 
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In our model, this component will be assumed to be 
of th? same magnitude for all the points c£ the cross 
section. 

Figure B-8 shows the resulting velocity diagrams for 


three particles located at radius r^, ^2 and rj. 



ItfV’Okiol direction 


I 



ot r»r, 

Fig. B-8. Velocity diagrams for particles at 
r^, X 2 and rj 
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It is ovident from these velocity diagrams that the 
fluid particles would move toward the fore end of the tank 
(the end toward which the screw seems to advance) follow- 
ing at different radii, different trajectories. 

Figure B-9 shows the velocity diagram at a radius 
r-0 (at the root of the helical screw, if the diameter of 
the screw shaft is disregarded) and at r-R^ (at the outer 
tip of the screw) . 



■ R 


s 


Fig. B-9. Velocity diagrams at r ■ 0 and r 
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The tangential coaponent of the velocity of the fluid 
with respect to the tank is given by 

V£/t cos a ■ V cot a (B8) 

for r ■ 0, while at r ■ is given by 

Vf/J stn e . ( . - Vf/, cos a (B») 

■ «R - V cot a (BIO) 

(since ? n a - cos 6.,, or at any radius by 

wr - V cot o (BID 

The distribution of these tangential velocity 
conponents, for the particles of fluid lyiug on a plane 
normal to the axis of rotation, is as shown in Figure 
B-10. 
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Fig. B-10. Distribution of the tangential 

velocity coaponents of the fluid. 


The velocity conponents for a particle of the fluid 
in cylindrical coordinates (r, 6, z) are shown in 
Figure B-11. The coordinate systen is fixed with respect 
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• u»r-Vcrt « 


^oxit 0? rototion 
^of fht tertw 


Fig. B-11. Velocity coaponentx for « particle 
of fluid in cylindrical coordinatet 
(r, 0, i) 


to the tank. The velocity field is thus described by 


V 


r 


• 0 


V 


e 


«r - V cot o 


V 


z 


- V 


(B12) 

CBIS) 


(B14) 
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Substitution of this velocity field into the Nevier- 
Stokes equations • vith the essuspiions that the fluid is 
inviscid and subjected to a zero g condition, yields 
the relations. 


2 

3p ■ p . Cdtr - V cotal 
3r r 


(BIS) 


iE " ® (B16) 

36 


and iE - 0 
3z 


(B17) 


FroB (B17), the pressure field,p, can be described 
as 

p ■ K ♦ ®) 


where 


K ■ constant 
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The &ixed derivfttives ^ P^(2r26) and ara 

identical, so that the given velocity field is a solution 
to the equations. 

Integrating the systen of equations results in: 


P 


^ Pw2 (r* - Rgi) 

Hg 

2 


(2 p 6) V cot u) (r-Rg) 


♦ (pV^ cot^ a) In 


r 


(Blf) 


where, Pg» is the pressure of the fluid at a radius Rg. 

ror a partially filled tank, since the pressure at 
the free surface of the liquid is unifom, the above 
equation indicates that the free surface of the liquid 
should be located at the sane radial distance along the 
tank. Purthemore, since the equation indicates that the 
pres ire Increases with r, considering that the pressure 
of tht gaseous phase is constant throughout the region it 


origw^ai: page w 
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occupies would limit the pressure distribution given by 
the above equation to the values of r for which R^r^Rg. 

Rs is the radius of the screw and Rg the radius at which 
the free surface is located. (Pressure variations within 
the gaseous phase would be negligible due to the order of 
magnitude of the density of the gaseous phase) . 

The above conclusions are illustrated in Figure B-12. 



Fig. B-12. Settling configuration 

If the mass of liquid is rotating as a rigid body, 
i. e. V - 0, then the pressure distribution is given by 


"i * 


P(0 


2 (j.2 


Rg^) 


2 


(B19) 
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In this case also, the free surface of the liquid 
stands as shown in Figure B-12. 

The above results points to an interesting design 
consideration in this model: the outlet for the liquid 

should be located at the periphery of the fore end of 
the tank. 



ORIGINAL PAGE IS 
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C. THERMODYHANIC CONSIDERATIONS 

From a thermodynamic point of view, the model is 
described by the system shown in Figure C-13. The liquid 
and gaseous portions within the tank are enclosed by a 
control surface. Work, in the form of paddle-wheel work, 
as well as heat, cross the boundary of this surface. 

Liquid is withdrawn from the control surface at a cons- 
tant rate while pressurant gas is added to maintain the 
pressure inside it constant. The pressurant gas is siq»plied 
at the temperature of the system. The latter is assmed 
constant during the process. 



Fig. C-13. Control volume 



-24- 


Application of the First Law of Theraodynanics to 
the systea results in: 

5W ♦ 6Q - da^ 


on as a rate equation as: 


dt dt 


dt 


* h. 
2gc ® 

♦ dap 

ad *h 

2gc 

• 

(B 0) 



n as: 



f ^ 

2gc • 

♦ 

dt 

V 1 

2«c 

■ m 


(Cl) 


~ (Bo) 
dt 


(C2) 


where, 


6W 

— ■ rate at which paddle-wheel work is supplied 
dt 


6Q ■ rate at which heat is added to the systea 
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■ rate of mss efflux (liquid) fr<ui the teak 
dt 


» rate of uass influx of pressurant gas into the 
dt 

tank 


vL. 

2gc 


kinetic energy of the fluid entering or leaving 
the control volune 


h ■ specific enthalpy of the mss of fluid leaving 
or entering the control voluM. 


^(Eo), rate of change in energy content within the 
dt 

control surface. 


(subscript a refers to conditions within the control 
volune. Subscript p refers to the state of the 
pressurant gas as it enters the control voluM and 
subscript e refers to the state of the liquid as it 
flows out of the control volune) . 
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Changes in gravitational potential energy are assumed 
zero. Local properties are assuaed constant during the 
process. 

The energy content within the control volume comprises 
internal energy U, rotational kinetic energy and trans- 
lational kinetic energy i. e., 

E 0 -fu ♦ Ep * ®k] 0 

SO that 

3t ■ 3t[" * ^ (C4) 

or 

!r • It f"!. • It W . • it fi ™ 

The internal energy within the control volume consists 
of the SUB of the internal energies of the liquid (U^) and 
gaseous (U^) phases in it, i. e., 

U - Uj ♦ Uj (C6) 
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The internal energy of the gaseous phase consists of 
the SUB of the internal energies of the propellant vapor 
and of the pressurant gas Up. 

Thus, 

U - ♦ Up (C7) 

- »1 ^1 ^ ®v ^ *p ^*p 

where, 

iDj^, and Bp ■ the Bass of the liquid, 
propellant vapor, and pressurant gas, 
respectively. 

and, 

Pi, , and Up ■ the specific internal 

energies of the liquid, propellant vapor, 
and pressurant gas, respectively. 
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Differentiating (C7), 

dU ■ d (Hj y^) ♦ d («y y^) ♦ d Mp) (C8) 

If the gaseous mixture is assumed saturated with 
propellant vapor, then at a given instant of time, its 
composition can be determined as follows: its pressure, 

p, is due to the partial pressures of the propellant vapor 
and pressurant gas. 

P * Py Pp (C9) 

also, 

y CVp) - (1-y) Vy (CIO) 

and 


y 


y (Wp) ♦ (i-y) Uy- 


(Cll) 


total pressure of the gaseous fixture 

partial pressure of the propellant vapor 
partial pressure of the pressurant gas 

specific voluae of the propellant vapor 

specific volume of the pressurant gas 

specific internal energy of the gaseous 
mixture 

specific internal energy of the propellant 
vapor 

spocific internal energy of the pressurant gas 

ratio of the mass of pressurant gas to the 
mass of the gaseous mixture, i. e., 
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From a table of properties for the propellant vapor, 
one obtains the following properties for the saturated 
vapor state, at the temperature, T, of the mixture: 


Py - (p sat.) Y 

Ccis) 

Vy ^^g^ T 

(C14) 

\ T 

(CIS) 


The partial pressure of the pressurant gas is then 

Pp • P ■ Pv (C16) 


Assuming that the pressurant gas behaves as an ideal 

gas. 


- RP T 


(C17) 


Now, 

y m ^ m 

l-y Vp (Rp T/pp) 


(Cll) 


Solving for y, 


1 ♦ &-I 

(Vg)^ Pp 


(C19) 


y is the ratio of the mass of pressurant gat to the mass 
of the mixture, so that (l*y)is the ratio of the mats of 
propellant vapor to the mass of the mixture. 

The gas constant of the mixture can be expressed at 

R ■ y Rp + (1-y) Ry CC20) 

where , 

Rp “ gas constant of the pressurant gas 
Ry * constant of the propellant vapor 

The mass of the gaseous mixture can be determined from 


m 


R T 


(C21) 
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where, 

■ volme of the ftseous fixture. 

The mass of propellant vapor is then 

- (1 - y) m (C22) 

and the mass of pressurant gas is 

«p ■ y m (C2S) 

Consider now the equilibrium composition of the gaseous 
mixture at a later instant of time. Assume that the 
efflux of liquid has been such that the volume of the 
gaseous mixture has increased to The composition of 

the gaseous mixture depends upon temperature and pressure. 
Since we are assuming them to be constant, the properties 
for the saturated vapor state at the temperature of the 
mixture remains constant. Following the procedure out- 
lined previously, we will conclude that y, and thus R, do 
not vary during the process. 
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The ntss of the gaseous mixture will depend upon the 
new vilune of the nixture, Vg2» 


*2 


P Vg2 

RT 


(C24) 


The mass of propellant vapor is now, 


m (1 ■ y) *2 

and that of the pressurant gas, 

• y »2 

The change in the mass of propellant \apor in the 
mixture is 


- (1-y) m2 * (1-y) m 
■ (1-y) (mj-m) 

[in] 


(C27) 
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snd that of the pressurant gas, 
- *P2 ‘ "P 


• y »2 ■ y « 
- y (««2 “ ") 



(C28) 


Thus, the mass influx of pressurant gas, as well as 
the change in the mass of the propellant vapor, are pro- 
portional to the change in volume of the gaseous mixture. 

The mass of liquid, m^ , that occupied the space 
v'Vg^ - Vg) is given by. 


(C29) 

''1 

where, 

“ specific volume of the liquid (taken as 
that of saturated liquid at the temperature 
T of the system) . 
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Of the quantity, anount Amy reaains in the 
gaseous phase, so that the mass of liquid which efflux 
from the system is 





\ 

I 

I 
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Using the relations » 

y • !!& (C2S) 

n 

■ f Vg2 - V») (C29) 

’'l 

and pV - m RT (C33) 

relation (C-3ii) , reduces to 

®efflux 
Recalling that 

''"p * y [ St ] 

The above can be sustituted in (C34) so that (C32) 
simplifies to 

-efflux - % [^' S] 

which shows that the influx of pressurant gas is propor- 
tional to the efflux of liquid. 
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Consider now the change In internal energy within 
the control volume: U ^ can be expressed as, 


Ug - Bj g, ♦ »P Up 

■ "l “1 * “v “v * "p j 

Multiplying and dividing the terms within the 
brackets by ( m^ ♦ Wp) , 

Ug - Bj Uj * 

“o ■ ">1 “l * 


(By ♦ Bp) 


( _!!£ ) Uv ♦ { — ^ 


Ay * Wp 


•v 


(»V 


* v[ 


(1-y) Vy ♦ y Pp 


] 


(C36) 



(C57) 


The term in the bracket now is y, the specific 
internal energy of the gaseous mixture, which is constant. 

Therefore, 



■ "l “I 


•f 


y 


(CSS) 
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Since depends upon teipereture and thus, is 
constant in our case, a change in the internal energy 
within the control surface is. 


AU^ ■ Amj ♦ y ♦ »p) 


+ y 


[*"v ♦ 

- y^ Anij ♦ y j Aaij (-I1-) ♦ A«j ( IL ) 1 

L Vp J 


All ■ Am, 


["i* 


w ( VI 


) * V Cil.) 




Substituting AjBj in terms of 


AU, 


<1 ^ P * V (]^) \ 


(CS9) 


■ efflux (C40) 


The terms within the brackets are constant, so that 
AU^ can be expressed as 


AUjj - Cj ^efflux 


where. 


C 


1 


yj ♦ y Vj 



(C41) 

(C42) 
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Consider now the change in rotational kinetic energy 
within the control voltune. The instantaneous rotational 
kinetic energy is deternined by sunaing the rotational 
kinetic energies for all the particles i within the 
control surface » i. e., 

2 2 

Ej. - E ”i ^fi ^i (C4S) 

2 gp 


For the fluid configuration shown in Figure C-14 one 
can determine the rotational kinetic energy of 





Fig. C-14. Fluid configuration for rotational 
kinetic energy considerations. 
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the particles for the section of the tank within radii, 
and R, by the relation 


(Ry R) 




where , 


Ry ■ radius of the screw shaft 

o)£ - angular velocity of the fluid particle 

If the contribution of the gaseous phase to the 
rotational kinetic energy is neglected and the liquid is 
assumed to remain during the process in the configuration 
shown in Figure B-12, then (C44) would represent the ro- 
tational kinetic energy lost due to the efflux of mass, 
considering that the tank was initially full. 
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The mass within the differential element of volume 
in Figure C-15, is given by 

dm “ p (Zwrdr) (L ♦ 2 t) (C45) 


where , 


L * length of the cylindrical portion of 
the tank 

l • variable length of the section at the 
hemispherical ends. 

The relationship between 1» R and r is given by 

9 

^ (C46) 


as shown in Figure C-IS. 
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Thus, 

dm ■ 2 irp L ♦ 2 v/r^^ - r‘ j r dr 
Dividing both sides by dt, 

2 TrpLr ♦ 4 ^Pt/r 2 - f2 1^ 

* Jdt 


dm 

dt 


(C47) 


(C48) 


gives the rate of change of the radius with respect to the 
rate of change of mass. It represents the rate of mass 
efflux of the liquid, “©ffiux* 


m 


efflux 


2TTpLr ♦ 4irpr 


dr 

a? 


(C49) 


Separating variables, 

"efflux “ ^2wpLr + 4iTpr y R^i - r^ j dr 
Integrating, considering constant, results in 

(csor 

which is the time it takes to reach a radius R, starting 


Tip 

^efflux 


at a radius R^. 
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Neglectins R^., 


I£_ 


m 


efflux 


CR,*-**)*'*) 


(CSl) 


From (C44) , the instantaneous rate of change of rota- 
tional kinetic energy is 


dF.r 

dt 


2 2 
r u)£ 


,2 , 2 

r u)£ 


dm 

3t 


efflux 


(C52) 


Substituting u»£ from (Bll) , results in 


dEr 

dt~ 


f2 r,i V cot g n2 ' 

C<*> r ' ®efflux 


(CSS) 


which is the quantity desired for the energy equation for 
the model under consideration. 


Regarding the rate of translational kinetic energy 
lost due to the efflux of mass, it can be determined from 

„2 . 


dEk 

dt 


dm 

2gc dt 


2*c 


'efflux 


(C54) 
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So far we have exaained the terms on the right hind 
side of the equation in (C2). ySn the left hand side of 
it, we have that 


jg>e 

dt 


"efflux 


(CSS) 


and is related to through 
"efflux " ^e ^e 


(CS6) 


where 

■ cross sectional area of the exit pipe 
p|^ • density of the liquid leaving the tank. 

I 

The term diUp/dt is related to through equation 

(CSS) , and is related to Vp through 

\ % Vp Pp 


(CS7) 
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where , 



cross sectional area of the pressurant 
gas inlet pipe 


Pp ■ density of the pressurant gis 


Substitution of the various relations derived into the 
energy equation (C2) would give the instantaneous relation- 
ship between power and heat transfer rate. 
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D. POWER REQUIREMENTS 

The assunption of an inviscid liquid in Section B 
wa» reasonable for the purpose of obtaining values for 
the velocities of the bulk of the liquid. On the other 
hand, a large fraction of the power requireaents are 
associated with viscous losses, and thus difficult to 
assess. 


Discussion of a method for predicting power require- 
ments for the model follows: If one considers thit the 

tank is initially full and no efflux of liquid occurs, 
while the latter is brought from rest tc the desired 
rotational speed, then the liquid would rotate approxi' 
mately as a rigid body and the energy requireaents for 
the liquid (disregarding that required by the helical 
screw itself) would correspond to the rotational kinetic 
energy of the liquid. It can be evaluated by integrating 


(C44) th R - R.. Thus. 


vp 


w2l R,^ 


4 

IT 


».*1 


CD-I) 




If the time required to achieve the desired rotational 
speed for the helical screw is t» then the average power 
requirement during this initial period would be 

? - (D-2) 
t 


Evidently, frictional and drag losses at the tip of the 
helical screw as well as inefficiencie: of the driving 
device will contribute additional power requirements. 

These can be better handled through tests and related 
through an initial settling efficiency as: 


ideal power requirement 
actual power requirement 


Once the screw has achieved the desired rotational 
speed, and the body of liquid rotates nearly as a rigid 
body, the fluid losses will be mainly due to viscous 
forces between the rotating body of liquid and the tank 
wall through the layer of liquid in the clearance space. 



An estimate of this loss can be obtained as follows: 
a cross-sectional view of the flow pattern at the tip of 
the screw is shown in Figure D-I6. The velocity of the 
fluid at the tip of the screw is given by 

U« - wRj - 2wNRg (D-4) 



Fig. D-16. Flow configuration between the tip 
of the screw and the tank waxi. 
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while that of the fluid at the tank wall It atsuaed 
zero. Due to viscous effects a velocit> distribution 
within the clearance space will be set up causing a 
shear stress and a corresponding drag force at the wall. 
The drag force can be estinated by applying the relations 
for flat plate boundary layers. 

The drag force D is given by 
a 

D-hJ^^Cx) dx (05) 

0 

where , 

b ■ length of the lay»>r in the radial direction 

X • length of the plate, in this case equal to 
the circumfevtjnce of the tank 


The ^hear stress at the wall, and corrsapondingly, 
the Method of analysis for D. will depend upon flow con- 
ditions. If. for example, laminar flow is assumed in the 
clearance space, the shear stress at the wall will be 
given by 


r 

'w 


,, ^ I 

^ I y ■ O 


(W) 
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where, 

y - perpendicular distance measured from 
the wall 

u ■ viscosity of the liquid 


Assuming a velocity distribution in the clearance 
space as» 


U - U* 




then. 




2M U* 


The drag force is then, 

2u 


2 wR^ b 


where. 


R. 


(D7) 


(D8) 


CD9) 


radius of the tank 
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Since the distance travelled per unit tim by the tip 
of the screw is 2nR^N, the power required to overcoae 
this frictional effect is 

p ■ 4 vv U«* b X 2ir Rj N (DIO) 

c 

2 

p ■ 8it R^ Rj M U« b N 
c 

Substituting for U*, results in, 

p - 16 ir kJ ub 
c 

Care should be used in applying this relation to the 
hewispherical ends. In this case a mean radius and the 
proper value of b should be used. 

In a siailar way, frictional losses of this type can 
be detemined in the case of turbulent flows. 

Now then, when the valve at the discharge end of the 
tank opens, allowing an efflux of liquid, a velocity 


(Dll) 


(D12) 
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of the fluid relative to the screw surface will develop. 
This velocity is related to V and the helix angle of the 
screw a by 


V - sin a (DIS) 

To determine the frictional losses in this case, one 
can assume the liquid to be moving along a helical channel 
of rectangular cross section formed by the walls of the 
helix, the tank wall and the screw shaft as shown in 
Figure D-17. 




Fig. D17. Channel-like conduit through 
the helical screw. 
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The standard procedure used for the detemination of 
losses in noncircular pipes can, thus, be applied. 

The cross sectional area, A^, available for the flow 

is 

\ - Rt* 

where , 

s ■ P cos a 
as shown in Figure D-18. 


P 
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The hydraulic dianater is given by 



4 X area 

wetted periaeter 


2 (R^-R) ♦ s 


(D14) 


Frictional losses in this case are deterained by the 
use of the Moddy Chart and standard procedures. 
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B. C0M4BNTS AND COMPITTATIONAL PROCEDURE 

It is expected that a simple model as the one 
previously discussed has shortcomings. In this particu- 
lar case, the shortcomings occur at both ends of the 
tank where boundary conditions are not satisfied by the 
velocity field chosen. This makes application of the 
model unsuitable for tanks with a small length to diame- 
ter ratio. It would be only applicable to large tanks, 
which probably comprise most of the storage vessels of 
practical use. 

At the fore end of the tank, one would expect that, 
although the fluid has been moving toward that end with 
a uniform velocity, V, not all of the liquid finds it way 
out of the tank. The portion of liquid remaining in the 
tank would recirculate locally or would try to find its 
way towards the aft end. This would result in large 
viscous losses. An approach that could reduce viscous 
losses in the model, with the added advantage that flow 
conditions would be closer to those assumed, would be to 
provide a return path for the liquid through the helical 
screw shaft as shown in Figure E-19. The return path 
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Flg. E-19. Inside of the screw shaft used as 
a return path for the liquid 
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would also provide an escape route for any gas trapped 
at the fore end during the initial settling of the 
fluid. 

The mass of liquid moving toward the fore end can 
be determined from 

m - Vp (El) 

where, 

- il - R^)Vp (E2) 

is the wetted portion of the screw’s cross sectional 

area. 

The rate of mass efflux “afflux 
tion K of m, i.e. , 

■efflux ■ * " 

(1-K) would represent the fraction of m which is 
recirculated. 

Thus , 

Vfflux ■ * "o 
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A sinuMiry of the coaputational procedure la pre- 
sented in whet follows: 

s) the desired rotational speed is detemined 
by the criteria set forth in Reference (3); 
the nininum rotational speed being jiven by 



b) the helical screw is started ;:md brcnght to 
the desired rotational speed. (A closed| full 
tank is assumed during this part of the pro- 
cess) . 

c) the discharge valve is opened and the desired 
liquid efflux is established. 

d) pressurant flow rate is determined from (CSS) . 

e) a small time interval is selected so that 
conditions can be considered quasi-steady 
during it. 

f) the radius R reached at the end of the time 
interval is determined by either (C50) or (C51). 
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ii 

g) the velue of V is deteniined froa (24) for an 
assuaed value of K. 

h) the pressure distribution is deterained using 
(B19). 

i) the change in internal energy within the 
control voluae is deterained froa (C41). 

)) the change in the rotational kinetic energy 
and the translational kinetic energy is 
determined using (CSS) and (C54) respectively. 

' k) frictional losses are determined by the 

procedures of section D. 

l) and Vp are evaluated from (C56) and (CS7) 
respectively. 

m) thermodynamic properties of the fluid are 
obtained from appropriate tables. 

n) all terms are substituted in the energy 

n equation and the instantaneous rate of heat 

transfer is determined. 
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o) the value of R for the next interval of tiae it 
deterained from (C50) or (CS1). 

p) the process is repeated for each interval 
of tiae during the transfer process. 
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P» SIMILARITY ANALYSIS AND MODEL TESTING 

Three types of similarity must be considered in 
settling systems. These are: geometric, kinematic 

and dynamic similarity. 

Two systems are considered geometrically similar 
when the ratios of corresponding diMnsions in one 
system to those of another have the same value. This 
implies that the systems, although of different sizes, 
have the same shape. 

Kinematic similarity exists between two systes» 
when, in addition to possessing geometric similarity, 
the ratios of velocities between corresponding points 
in each system have the same value. 

Dynamic similarity exists between two systems when, 
in addition to being geometrically and kinematically 
similar, the ratio of forces between corresponding points 
have the same value. 
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The diaensionless groups pertinent to the settling 
systen being considered are: 

2 


N 


Re 


P nd. 


(FI) 


2n.3 


N - P 


(F2) 


N 


Fr 


D*n2 


g 


(F3) 


^ • lo;- 


(F4) 


N. 




(F5) 


B 

as well as the geometric ratios X • 'Dg* 
Og/^ Ds/L(y Dg/Dfy nnd Ds/1^» 

where , 

Nge * Reynolds nuaber 

■ Weber nunber 

WC 

Npr ■ Froude number 

J ■ advance ratio 
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N^i - power nuiri>er 

N - rotational speed of the helical screw 

D, • diameter of the helical screw 

Dt " dianeter of the tank 

c ■ clearance between the tip of the screw 
and the tank wall 

p ■ density of the liquid 

V • viscosity of the liquid 

P* - power supplied to the systen 

P ■ pitch of the helical screw 

L( • length of the tank 

a - surface tension of the liquid 

" dianeter of the helical screw shaft 

■ diameter of the liquid discharge line 

Dp ■ diameter of the pressurant line 

g • gravitational acceleration 

axial velocity of the fluid (a function 
"efflux^ 


V 
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The Weber end Froude nuaber can be coiid>ined to for* 
the Bond nuiber, i.e,» 

Bo • ^ (P^) 

^ a 

Equality of all the groups assures similarity between 
systems of different size. The last seven dimensionless 
groups of the above list represent the groups necessary 
for geometric similarity which requires that all corre- 
sponding dimensions between a model and its prototype 

bear the same ratio to each other. The reference dimen- 
sion used for the above group is the helical screw diaiw- 

ter. 


Confining attention to geoMtrically similar systems 
would require equality of the following groups to insure 
dynamic and kinematic similarity: 


N 


Re 


N. 


we 


N 


Bo 


N, 


J 


An experimentally derived function correlating the 
data for such a system would be of the form 


f( 


P ND. 


2 3 

P 


P 

a 


pnSdjS '*“7 


) - 0 (F7) 
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Consider the subscript n to represent the value of 
a property for a model while subscript p represents the 
value of the same property for its prototype. 


Model testing requires that the values of correspond- 
ing dimensionless groups between model and prototype be 
the same. Test conditions that must be satisfied to 
insure similarity of solutions will be analyzed in what 
follows. To begin, consider the Froude number. To sat- 
isfy similarity conditions requires that the ratio of the 
Froude number of the prototype to the Froude number of 
the model satisfy 


(Np,)p 


1 


(F«) 


or 




sm 



(^ ) - 1 


(F9) 


Solving for Nm/Np yields 



(FIO) 


Considering 



as 


10^ and ®sp^®sm as 30, then 


^ - (30)^ (10^) - 9 X 10® 


(F11) 
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If the prototype were to turn at 6 rpa^ the above 
relation Means that testing at one g would require a 
rotational speed of 

- S.4 X 10® rpii CF12) 

which is an excessively high value. 

However » if one assuMs that the naxiMua rotational 
speed of the screw is 84 rpn, then, the Maxima g level at 
which testing would be satisfactory would be. 



(F13) 


- 1.6 X 10 

Preservation of similarity also requires that the 
ratio of Reynolds number also be 1. Thus, 



1 (F14) 
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so that 



(FIS) 


If one considers water at 1 atai, 300K (^«540R) as 
a referee fluid, and cryogenic liquid hydrogen saturated 
at 24. 7K (4SR) as the prototype fluid, then 



(F16) 


(FI 7) 


Substituting these values in the above expressions 
results in 


which is much larger than the assuised ratio of 14. This 
■eans that a different referee fluid under the proper 
therwodynanic state must be sought. 


Considering now the Weber number, again, 


TO 




so that. 


(F19) 



(P20) 


For the water-hydrogen combination considered 
previously, 

r ■ 59 (F21) 

P 

and. In this case. 


Nm _ 

IT ■ tTHF’ 

P 

- 320 


(F22) 


which Is larger than the assumed value of 14, pointing 
again to the need for a proper referee fluid. 

The power number is useful in predicting power 
requirements foi the prototype once the experimental values 
for the model are obtained. 
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Thus, 



ky[‘t.y 


s 


•1 


(F2S) 


(F24) 


The ratio of J^/J. is useful in determining proper 

P nl 

ratios of flow rates for the liquid discharge since 
"efflux proportional to V. 
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GONCLUSION 

An initial attempt to provide a nodel for engineer* 

ing analysis of a screw settling device used as a iMsans 
of settling the liquid-vapor phases of a fluid within 
a tank for subsequent transfer to another vessel, has 
been presented in this work. Actual tests will indicate 
the suitability of the nodel, or provide the basis for 
necessary Modifications. 

Similarity analysis indicated that: 

a) it is impractical to carry out tests 
at one g. 

b) an adequate referee fluid should be 
sought for tests with a scaled down 
model of the system. 
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